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Major Drivers of Current Work 






1 To provide computational tools as an economical option for developing 
future space transportation systems (i.e. RLV subsystems development) 

Impact on component design Rapid turn-around of high-fidelity analysis 

Increase durability/safety => Accurate quantification of flow 

(i.e. prediction of flow-induced vibration) 

Impact on system performance => More complete systems analysis 

using high-fidelity tools 

• Target 

Turbo-pump component analysis Entire sub-systems simulation 
Computing requirement)^1^g^\ 

=^>The goal is to achieve 1000 times speed up over what was possible in 1992 


■y> 
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Current Challenges 


• Challenges where improvements are needed 

- Time-integration scheme, convergence 

- Moving grid system, zonal connectivity 

- Parallel coding and scalability 

• As the computing resources changed to parallel and distributed 
platforms, computer science aspects become important. 

- Scalability (algorithmic <& implementation) 

- Portability, transparent coding, etc. 

• Computing resources 

- "Crid" computing will provide new computing resources for 
problem solving environment 

- High-fidelity flow analysis is likely to be performed using "super 
node" which is largely based on parallel architecture 
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Time Accurate Formulation 


Time-integration scheme 

Artificial Compressibility Formulation 

■ Introduce a pseudo-time level and artificial compressibility 
* Iterate the equations in pseudo-time for each time step until 
incompressibility condition is satisfied. 


Pressure Projection Method 

• Solve auxiliary velocity field first, then enforce 
incompressibility condition by solving a Poisson equation 
for pressure. 
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INS3D - Incompressible N-S Solver 


Ms Smsrti tester^ 


** Parallel version is based on INS3D- UP 
•MPI and ML P parallel versions 

• Structured, overset grid orientation 

• Moving grid capability 

• Based on method of artificial compressibility 

• Both steady-state and time-accurate formulations 

• 3 rd and 5 +h -order flux difference splitting for convective terms 

• Central differencing for viscous terms 

• One- and two-equations turbulence models 

• Several linear solvers : CMRES, GS line-relaxation, LU-SGS, 

GS point relaxation, ILU(0)..,.. 

•HISTORY 

** 1982-1987 Original version of INS30 - Kwak, Chang 
** 1988-1999 Three different versions were devoped : 
INS3D-UP / Rogers, Kiris, Kwak 
INS3D-LU / Yoon, Kwak 
INS3D-FS / Rosenfeld, Kiris, Kwak 
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Tim® (see) p®r iteration 
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RL V 2 nd Gen Turbopump 


VELoary magnitude 
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- 34.3 Million Points 
-Two and half impeller 
rotations are completed. 

♦One complete rotation 
requires less then 3.5 days 
by using 128 CPUs on SCI 
Origin 3000. When 512 
CPUs are utilized one 
rotation can be completed 
less then 1.5 days. In 1999, 
one impeller rotation would 
take 42 days by using 32 
CPUs on SCI Origin 2000 
platform. 


RLV 2 nd Gen Turbopump 
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L21CPtJsG3K<40QMhl) 


MLWckihe optimized 
12J CPUs p3K(400Mffe) 


May-2000 


! 1.3 days 

MLiP/Cache optimized 
480 QKJt 031C(400MHzi 

* i i i — i — l 

May-2001 Sep-' 
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Scripting Capability 




SCRIPT DEVELOPMENT FOR TURBOPUMP SIMULATIONS 
Motivation 

Significant user’s effort needed 
in complex process from geometry ’SV 
to flow solver , mp e 

Objective 

Develop script system to 

- generate grids 

- create domain connectivity Input 

- create flow solver Input 

for different components diffuser-^ ^ . 

automatically turSojMibip Scfifts Spi^ 


Approach 


(grids I domain connectivity Input I |ffow aohrar input] 


Develop one script for each component with ring Interface 
between components => easy plug In for different designs 
and combinations of components 
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Scripting Capability 


SCRIPT GENERATION 

Disadvantages 

> Require expertise to build scripts the first time 
Advantages 

> Allow rapid re-run of entire process 

> Easy to do grid refinement and parameter studies 

> Easy to try different grlddlng strategies 

> Documentation of grlddlng procedure 
Tel scripting language 

> Works on UNIX, LINUX and WINDOWS 

> Integer and floating point arithmetic capability 

> Modular procedure calls 

> Easy to add GUI later If needed 
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Scripting Capability 



INPUT* AND OUTPUT 

Input 

> profile curves for hub and shroud In PLOT3P format 
(rotated by script to form surface of revolutlo n) 

> blade and tip surfaces in PLOT3D format 


> Parameters that cart be changed 

- number of blades and sections 

« global surface grid spacing As (on smooth regions) 
local surface grid spacing, sbme Independent and some 
expressed as multiples of as (feadl ng/tra III ng edges, etc.) 

- normal Wa II grid spacing (viscous, Wall function) 

- marching distance 

- grid stretching ratio 


Output 

> overset surface and volume grids for hub, shroud, blades 

> object X-rays for hole cutters using DCF 

*> domain connectivity namelist Input for 0VERFI.OW-P 


31 






Scripting Capability 
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IMPELLER 

(M sections, N different blades In each section, fljp clearance) 








^ V & 
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Manual Script (fine) Script (coarse) 
No. of pts (million) 19.2 l]5.2 8,8 

User time* ~ 2 weeks 31 9 sec. 234 sec. 

c from geometry del. to DCF fnputwfth $GI R12k 300MHz CPU) 


Scripting Capability 


DIFFUSER (N repeated blades, no tip clearance) 


[NASA 











Manual 

Script (fine) 

Script (coarse) 

No. of pts (million) 

8.6 

1,6 

User tfme * 1 day 

{* from geometry def. tp DCF input Witt) 

37 sec. 

SGI R12k 300MHz CPU) 

22 sec. 
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STATIC/DYNAMIC STRESS ANALYSIS 
FOR TURBOPUMP SUB-SYSTEMS 


NEW CONDITIONS 

, INTERFACE 
ZIPPER/GRID ” 


CFD GRID T 

PRESSURE I 

TEMPERATURE FEM GRID 


STRUCTURAL 

LOADS 


FLUIDS 

INS3D 


STRUCTURES 

NASTRAN/ANSYS 


STRESS 

ANALYSIS 


Summary 
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•Unsteady flow simulations for RLV 2 nd Sen baseline turbopump for three 
impeller rotations are completed by using 34.3 Million grid points model. 

•MPI/OpenMP hybrid parallelism and MLP shared memory parallelism has 
been implemented in INS3D, and benchmarked. 

• For RLV turbopump simulations more than 30 times speed-up has been 
obtained. 

•Moving boundary capability is obtained by using DCF module. 

• Scripting capability from CAD geometry to solution is developed. 

• Fluid/Structure coupling is initiated. 
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